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Abstract 

We present a tight-binding study of donor impurities in Si, demonstrating the adequacy of this 
approach for this problem by comparison with effective mass theory and experimental results. We 
consider the response of the system to an applied electric field: donors near a barrier material and in 
the presence of an uniform electric field may undergo two different ionization regimes according to 
the distance of the impurity to the Si/barrier interface. We show that for impurities ~ 5 nm below 
the barrier, adiabatic ionization is possible within switching times of the order of one picosecond, 
while for impurities ^ 10 nm or more below the barrier, no adiabatic ionization may be carried out 
by an external uniform electric field. Our results are discussed in connection with proposed Si:P 
quantum computer architectures. 

PACS numbers: 71.15.Ap, 71.55.-i, 71.55.Cn, 03.67.Lx 
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I. INTRODUCTION 



Simple donors in Si have recently become the subject of renewed interest due to proposals 
of quantum computer architectures in which P donors in Si play the role of qubits.i^i^ Logic 
operations in such architectures involve the response of the bound electron wavefunctions to 
voltages applied to a combination of metal gates separated by a barrier material (e.g. Si02) 
from the Si host. The so-called A-gate, placed above each donor site, pulls the electron 
wavefunction away from the donor, aiming at partial reduction^ or total cancelation^ of 
the electron-nuclear contact coupling in architectures where the qubits are the ^^P nuclear 
spins. In a related proposal based on the donor electron spins as qubits,— the gates drive 
the electron wavefunction into regions of different (/-factors, allowing the exchange coupling 
between neighboring electrons to be tuned. Ideally, electric-field control over the donor 
electron wavefunction requires all operations to be performed in the adiabatic regime^, which 
sets a lower bound for the time scales involved in such processes. 

Recent studies have demonstrated that the tight-binding (TB) approach, traditionally 
adopted for deep levels,— provides a valid description for intermediate^ and shallow levels^ 
in semiconductors. Impurity states are calculated from a sequence of supercell sizes and a 
finite-size analysis which provides extrapolation to the bulk limit. Also, electric-field effects 
may be easily incorporated within the TB scheme^, allowing estimates of switching times in 
electric-field-tunable devicesA*?-^. In this work we present a TB description for donors in Si, 
aiming at a physical description of the relevant properties involved in the A-gate operations 
mentioned above. 

Donors in Si have been extensively and successfully investigated within the effective mass 
theory (EMT),— thus providing a preliminary test for the TB approach by comparison of 
wavefunctions predicted by the two formalisms. This comparison is presented in the next 
section. In Sec. IIIII we explore a simplified model of the A-gate operations in the Kane 
quantum computer proposal^ by considering the Si:P system under an uniform electric field 
and near a barrier. In Sec. II VI we discuss operation times and restrictions imposed by the 
donor positioning with respect to the Si/barrier interface in connection with the adiabaticity 
of the A-gate operations. Our summary and conclusions are presented in Sec. El 
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II. TB DESCRIPTION FOR DONORS IN SILICON 



A. 



Formalism 



The TB Hamiltonian for the impurity problem is written aa^ 



(1) 




where i and j label the atomic sites, /i and v denote the atomic orbitals and rj is the distance 
of the site i to the impurity site. The matrix elements h'^^ define all the on-site energies and 
first and second neighbors hoppings for the bulk material. The donor impurity potential 
U{r-i) is described by a screened Coulomb potential (e = 12.1 for Si) 



At the impurity site (r^ = 0), the perturbation potential is assigned the value — f/o, a param- 
eter describing central cell effects characteristic of the substitutional species. In the present 
calculations, f/o was kept as an adjustable parameter (previous estimates for this parameter- 
are of the order of one to a few eV). We adopt here the sp^s* TB parametrization for Si pro- 
posed by Klimeck et al,^ which includes first and second neighbors interactions. Inclusion of 
hopping matrix elements up to second neighbors provides a good description of the effective 
masses at the conduction band minima. This parametrization gives the fc-space positions of 
the six band minima at the six equivalent points along the A lines, at Amin = 0.75(27r/a), 
where a = 5.431A is the conventional cubic lattice parameter for Si. We do not include 
spin-orbit corrections in our calculations. 

The eigenstates of H are determined for a system where a single impurity is placed in 
a cubic supercell containing = atoms arranged in the diamond structure, where L 
is the length of the supercell edge in units of a. The supercells are subject to periodic 
boundary conditions, and full numerical diagonalization can be performed for L < 6. Much 
larger supercella^ (up to 10^ atoms) may be treated within a variational scheme^^ where the 
ground state wavefunction and binding energy El for a donor level is obtained by minimizing 
the expectation value of (\E' | (if — £re/)^| For the donor ground state, Sref is a reference 
energy chosen well within the gap, but nearest to the conduction band minimum, and excited 
states are obtained by tuning e^e/ towards the conduction band edge. Finite size scaling 




(2) 
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allows extrapolation to the bulk limit {L 



oo) according to the ansat^^ 



(3) 



where Ei, is the binding energy for a single donor in the bulk. 

The eigenf unctions of (Q) written in the basis of atomic orbitals |0i/(r — Rj)) are given 
by |\l/TB(r)) = J2iu^i>^\^'^(^ " -^i)) where the expansion coefficients ai^, give the probability 
amplitude of finding the electron in the orbital u localized at Rj. We do not include explicit 
expressions for the atomic orbitals; the overall charge distribution is conveniently described 
through the TB envelope function squared,— 



B. Donor ground state 

In the proposed TB model, the only free parameter is related to the on-site value for the 
impurity potential, Uq. In Fig. ^a) we present the converged (L —>■ oo) binding energy of 
the lowest donor state as a function of Uq. We also characterize the donor ground state by 
its orbital averaged spectral weight^ at A^m 



where N is the number of atomic sites in the supercell, and the first summation is over the 
six equivalent at the conduction band minima. This quantity is plotted in Fig. ^b) as a 
function of Uq. 

We determine the value of Uq so that the binding energy of the donor results to be in good 
agreement with the experimental value which, for P in Si, is Eb = 45.6 meV. As indicated 
in Fig. ^ Uq = Up = 1.48 eV gives the correct binding energy for the P donors in Si. This 
value for Uo is used in the calculations below. 

C. Comparison with EMT 

EMT exploits the duality between real and reciprocal space, where delocalization in real 
space leads to localization in fc-space, e.g. for shallow donors around the fc-vector at the 
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minimum of the conduction band. Within EMT in its simplest formulation,— the ground 
state for donors in Si is six-fold degenerate, due to the six-fold degeneracy of the Si conduc- 
tion band. Valley-orbit interactionsi^ lead to a non-degenerate ground state wavefunction 
of y4i symmetry,— 

^W = 7^E^M(r)0M(r), (6) 

V fj.=i 

where ^/^(r) = M^(r)e*''f '"^ are the pertinent Bloch wave functions, and the envelope functions 
given by (e.g. for /i = z)^ 

The effective Bohr radii for Si from a variational calculation are a = 2.51 nm and b = 
1.44 nmM In Fig. |21we present the TB envelope function squared calculated from (jU along 
three symmetry directions with the corresponding EMT results obtained from (jHl), where 
the periodic part of the Bloch functions have not been explicitly included, consistent with 
not explicitly including the atomic orbitals in the TB description. Note that the oscillatory 
behavior coming from the interference among the plane-wave part of the six 0^ is well 
captured by the TB envelope function. 

The good agreement between TB and EMT is limited to distances from the impurity site 
larger than a few lattice parameters (~ 1 nm). Closer to the impurity, particularly at the 
impurity site, the TB results become much larger than the EMT prediction, in qualitative 
agreement with experiment.— This reflects central cell effects, not included in the EMT 
expressions © and (|7j). In the central cell region, the discrepancy between TB and EMT 
wavefunctions is significantly larger than those reported for donors in GaAs,- a result that 
could have been anticipated from the spectral weight given in Fig. ^b). EMT rests on the 
assumption that the impurity eigenstate is highly localized in fc-space, so that only Bloch 
states near the conduction band minima enter in the expansion, as implied in Eq. (jH)). This is 
the case for GaAs,— where for a range of values of Uq {Uq < 1.8 eV) we find W{T) essentially 
equal to one, in agreement with the EMT assumption. In Si, even small values of Uq yield 
spectral weights at A^m well below one. For Uq = Up in particular, W{Amin) — 0.3. 

We remark that the sharp shallow-to-deep transition obtained for GaAs in Ref. iSi, with 
kinks in the curves of Ef, and W versus Uq, is not reproduced here (see Fig.[H). We attribute 
this to the lack of a strictly shallow region, with the spectral weight of the donor state 
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concentrated in one or a few fc-points. Therefore, while the binding energy of shallow 
donors in GaAs is essentially constant, independent of the species (~ 6 meV for C, Si and 
Ge, in excellent agreement with the EMT estimate), in Si it varies according to the donor 
species (45 meV for P, 53 meV for As and 42 meV for Sb, to be compared with the EMT 
estimate of 30 meV). It is interesting to note in Fig.^a) that, as the impurity level becomes 
shallower by decreasing Uq, Ef, approaches the EMT estimate for the binding energy.— 

III. DONORS UNDER AN UNIFORM ELECTRIC FIELD 

The formalism presented in SecOis easily extended to include an uniform electric field in 
the system. Assuming a constant field E applied along the [001] direction, it is incorporated 
in the TB formalism by modifying the on-site energies in Eq. (P) as follows:^^*^^ 

h'^riE) = /.-(O) - \e\Ez.. (8) 

Periodic boundary conditions lead to a discontinuity in the potential at the supercell bound- 
ary Zi = Zb, where Zb is half of the supercell length along [001] or, equivalently, the distance 
from the impurity to the Si/barrier interface. The potential discontinuity, Vb = 2\e\EZB, 
actually has a physical meaning in the present study: It models the potential due to the 
barrier material layer above the Si host^ (see inset in Fig. E)). 

A description for the A-gate operations may be inferred from the behavior of the TB 
envelope function squared at the impurity site under applied field E, normalized to the 
zero- field value: 

A/Ao = \^Uo)\V\^%Fm' (9) 

The notation here indicates that this ratio should follow a behavior similar to that for the 
hyperfine coupling constants between the donor nucleus and electron with (A) and without 
(Aq) external field. Since the hyperfine interaction A is proportional to |\I'(0)p, and we 
are using here the envelope rather than the full TB eigenfunctions, this equivalence is not 
rigorous. The ratio in © is plotted in Fig. Efa) for three values of the impurity depth 
with respect to the Si/barrier interface. Calculations for ^^=10. 86 nm were performed 
with cubic supercells {L = 40), while for Zb= 5.43 and 21.72 nm tetragonal supercells 
with Lx = Ly = 40 and = 20 and 80 respectively were used. At small field values we 
obtain a quadratic decay of A/Aq with E, in agreement with the perturbation theory results 
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for the hydrogen atom.— At large enough fields, |\l/f^(0)p becomes vanishingly small, and 
the transition between the two regimes is qualitatively diff'erent according to Zb'- For the 
largest values of Zb we get an abrupt transition at a critical field E'c, while smaller Zb (e.g. 
Zb = 5.43 nm) lead to a smooth decay, similar to the one 

case, we define E^. as the field for which the curve A/Aq vs E has an infiection point, where 
A/Aq ~ 0.5, thus £'c(5.43nm) = 130 kV/cm. We find that the decrease of Ec with Zb 
follows a simple rule Ec oc I/Zb-, as given by the solid line in Fig. |2fb). 

In order to analyze the different regimes illustrated in Fig. Efa), we study the overall 
behavior of the envelope squared profile along the z-axis, 



where the first summation is over the two fee sublattices, with corresponding to the 
atomic sites in sublattice s, thus z labels each monolayer in the diamond structure, and 
p{z) quantifies the z-projected charge distribution for the electron states under applied field 
E. Fig. m gives p{z) for the electron ground state and also for the first excited state with 
Zb = 10.86 nm as the applied field increases. Up to fields very close to E^. (~ 53 kV/cm), 
the ground state distribution retains essentially the bound donor character, with the elec- 
tronic charge accumulating predominantly around the impurity [z = 0). For E > Ec we 
observe an abrupt charge transfer towards the barrier, with some residual charge remaining 
at the impurity site. The first excited state displays a complementary behavior, with charge 
transfer from the barrier into the impurity region as E increases. The binding energies (en- 
ergy eigenvalues relative to the bottom of the conduction band) are calculated here taking 
into account the dependence of the conduction band edge under applied field. The binding 
energies of the two lowest electron states are given in Fig. El^a). Note that they cross at Ec- 
The binding energies of the two lowest eigenstates for Zb = 5.43 nm are presented in 
Fig. Efb). They do not cross, but rather display an anticrossing behavior, confirmed by the 
corresponding doubled-peaked charge distributions in Fig. IHl with wavefunctions extending 
over the attractive wells of the impurity and of the electric field potential. This is consistent 
of eigenstates which are superpositions of bound states in each potential well. Note that for 
E = Ecin Fig.infc), the two states have essentially the same charge distribution, as expected 
at the anticrossing point. The anticrossing in Fig. Efb) is such that for E < Ec the lines 
giving the two states are essentially parallel, converging asymptotically at zero field to the 
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binding energies 45.6 meV, for the Ai ground state, and 32.4 meV for the first excited state. 
This is very close to the experimental binding energy of the excited E (32.6 meV) and T2 
(33.9 meV) states, which can not be individually resolved within our variational scheme.— 
Note that this was independently obtained with the same value of the parameter Uq, chosen 
to fit the Ai state binding energy alone. Near and above a typical 2-level anticrossing 
behavior is obtained, with the excited state eventually merging into the conduction band at 
E = 150 kV/cm. 

The above results may be understood within a simple picture of the electron in a double 
well potential, the first well being most attractive at the impurity site, V{Ti = 0) = —Uq, 
and the second well at the barrier interface, V{z = Zb) = — Vb/2 = —\e\EZB neglecting 
the Coulomb potential contribution (0) at the interface. An internal barrier separates the 
two wells and, for a fixed E, this internal barrier height and width increase with Zb- Deep 
donor positioning leads to a weaker coupling between the states localized at each well, even 
close to level degeneracy, resulting the level crossing behavior illustrated in Figl^a). For 
donor positioning closer to the interface the internal barrier gets weaker, enhancing the 
coupling between levels localized in each well and leading to wavefunction superposition and 
to the anticrossing behavior illustrated in Figl^l^b). The scaling of E^. with 1/Zb niay also 
be understood assuming that the critical field corresponds to the crossing of the ground 
state energies of two wells: The Coulomb well and an approximately triangular well at the 
barrier. Since the relative depths of the wells increases with EZb, and assuming that the 
ground states energies are fixed with respect to each well's depth, leads to the Ec oc 1/Zb 
behavior. 

IV. ADIABATIC PROCESSES DRIVEN BY AN UNIFORM ELECTRIC FIELD 

Coherent manipulation of electrons by the A-gates requires that the switching time be- 
tween different electron states be slow enough to guarantee adiabaticity of the process. 
Instantaneous eigenstates of H{t) may thus be defined at any time t. In the present case, 
we assume a linear increase of the external field from to a maximum value -Emax so that 
H{t) = H{0) — |e|i?max-2^, with < t < T, where T is the total switching time. A lower 
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bound for T is obtained from the adiabatic theorem,i^^ following RefllOt 

T„ = (11) 

9 mm 

where Qmin is the minimum gap between the two lowest electron states. In the anticrossing 
case illustrated in FiglSfb), we get g^i^ = 9.8 meV. Assuming that a totally ionized state 
is required as the final state, we take -Emax = 180 kV/cm, leading to Ta ~ 0.5 ps. This 
is a perfectly acceptable time for the operation of A-gates in spin-based Si QC, given the 
relatively long electron spin coherence times (of the order of a few ms) in Si.— 

As the impurity distance from the barrier increases, one eventually reaches the crossing 
regime, when Qmm —>■ 0, meaning that — > oo and no adiabatic ionization is possible. 
Ionization would still occur for E > Ec, but as a stochastic decay process from the first 
excited state. From Fig. IHl^a) we see that the A-gate might be used to partially reduce the 
contact interaction, in the case of Zb = 10.86 nm to about 20% of its value at zero field. For 
larger Zb the range for adiabatic variation in A/Aq is even smaller. Therefore Zb ~ 5 nm 
seems to be a favorable positioning for the donors, since it allows adiabatic reduction of 
A/Aq to any desired final value, with this ratio varying smoothly from one (at E = 0) to 
zero (for E = ^^ax ~ 2Ec). 



V. SUMMARY AND CONCLUSIONS 



We have presented a TB study of donor levels in Si. The reliability of the TB approach 
for the present study was verified by comparison of the TB and EMT envelope functions 
as well as by the value predicted for the Ai — {E,T2} energy splitting in agreement with 
experiment within our numerical accuracy. Previous TB studies of intermediate and shallow 
impurity levels in semiconductors^!^ dealt with materials with band extrema at k = 0, and 
the present results show that the oscillatory behavior of the wave function due to interference 
effects in the plane-wave part of the Bloch wavefunctions, typical of degenerate band extrema 
at k 7^ 0, is well is well captured by the TB approach. 

In the presence of an increasing uniform electric field, the donor states respond in different 
ways according to the donor depth Zb below the Si/barrier interface. For deeply positioned 
donors, i.e, for Zb >> ajb, where a and b are the Bohr radii for P in Si, abrupt ionization 
occurs at a critical field E^, while for Zb greater but of magnitude comparable to the Bohr 
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radii, a smooth electronic charge transfer from the donor site towards the barrier interface is 
obtained, eventually leading to complete ionization. The different regimes were identified in 
three ways: (i) From the decrease in electronic charge at the donor nucleus [Fig. |21^a)]. This 
behavior implies an analogous dependence of the electron-nucleus hyperfine coupling con- 
stant A as a function of the increasing external field, (ii) From charge distributions (Figs. E] 
andini), where the superposition of donor-like and barrier-like bound states is inhibited for 
deeply positioned donors, (iii) From the behavior of the binding energies of the two lowest 
electron states as the applied field increases (Fig. E)), changing from a level-crossing into an 
anticrossing regime as Zb decreases. The donor excited states in the S-like manifold also 
play a role in the anticrossing regime, as illustrated for < E'c in Fig-IMt')- 

The minimum gap g^am in the anticrossing regime is a key ingredient determining the 
possibility of an adiabatic evolution of the electron state under the action of the A-gates. 
Given that the product E^Zb is approximately constant [see fit Fig. El^b)], the adiabatic time 
Ta in (|TT| is expected to depend very weakly in the product E^s.xZb, assuming one aims at 
complete ionization.i^ Therefore Ta should not depend explicitly on Zb, but only implicitly 
through We have shown that for Zb ~ 5 nm, i.e., about twice the largest Bohr 

radius a in Eq. ((7j), electric field switching times smaller than 1 ps may be reached, which is 
a favorable operation time given the long electronic spin coherence times in Si. If one aims 
at a final state where only partial reduction of the electronic charge at the nucleus occursi*^, 
values of Zb of this order of magnitude are still the most convenient, since any final value 
of the nuclear charge may be attained. 

The Bloch phases interference behavior in the donor wavefunctions has been previously 
shown to lead to oscillatory behavior of the exchange coupling between two donors, affect- 
ing the two-qubit operations in exchange-based architectures in Si. We remark that such 
oscillations are well captured in the TB wavefunctions, and that the present study demon- 
strates that electric field control over single donor wavefunctions, such as proposed in A-gate 
operations,'^*^ do not present additional complications due to the Si band structure. The 
only critical parameter is the donor positioning below the Si/barrier interface, which should 
be chosen and controlled according to the physical criteria presented here. 
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FIG. 1: (a) Binding energy of the ground impurity state as a function of the on-site perturbation 
strength Uq, obtained from the L — *■ oo extrapolation ansatz. The dotted hne indicates the value 
Uq = Up that reproduces the experimental Si:P Ai state binding energy, (b) Calculated spectral 
weight at the conduction band edge for the ground state. Notice that as the perturbation Uq 
becomes weaker, Ej, approaches the EMT binding energy, while W does not approach 1. 
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FIG. 2: The dots give the TB envelope function squared for the lowest impurity state along three 
high-symmetry directions. The lines are the corresponding effective mass IV'I^ results. Note that 
the TB approach captures the oscillations of the EMT wave function in the asymptotic region. 
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FIG. 3: (a) TB envelope function squared at the impurity site under applied field E, normalized to 
the zero-field value, for the indicated values of the impurity-Si/barrier interface distance Zb- (b) 
Dependence of the critical field Ec on Zb- The solid line is a best fit of the form Ec oc 1/Zb- The 
inset gives a schematic representation of the perturbation potential added to the bulk Si hamiltonian 
due to the impurity at R = and to a uniform electric-field in the negative z direction. 
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FIG. 4: Tight binding envelope function squared projected along z for Zb = 10.86 nm and the 
indicated values of the field E applied in the negative z direction. The solid (dashed) line gives 
the donor ground (1st excited) state. Note in (b) and (c) the exchange among the p{z) for the 
lowest energy states [ground\ [excited\ which occurs over a narrow rage of electric field increase, 
a signature of the crossing behavior in Fig. [SJ^a). 
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FIG. 5: Calculated binding energies versus electric field intensity of the two lowest donor electron 
states, (a) For Zb = 10.86 nm the energies reveal a crossing regime, (b) Anticrossing of the two 
lowest electron states for Zb = 5.43 nm. The open symbols correspond the zero field calculated 
values: 45.6 meV and 32.4 meV. 
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FIG. 6: TB envelope function squared projected along z for Zb = 5.43 nm and the indicated values 
of the applied field E. The solid (dashed) line gives the ground (1st excited) state. At the critical 
field in (c) the two states have similar charge distributions, typical of a superposition of states 
localized in each well and a signature of the anticrossing behavior in Fig. [Sfb). 
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